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Abstract 



We study how the entropic uncertainty relation for position and momen- 
tum conjugate variables is minimized in the subspace of one-dimensional an- 
tisymmetric wave functions. Based partially on numerical evidence and par- 
tially on analytical results, a conjecture is presented for the sharp bound and 
for the minimizers. Conjectures are also presented for the corresponding sharp 

Hausdorff-Young inequality. 
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Let i/j be in a square integrable function in R n , to represent the wave function of a 
quantum-mechanical particle, and let p be its normalized probability density, to wit, p(x) = 
|í/'(x)| 2 /||?/'||2, where \\ip\\ p denotes the p-norm (/ \íjj(x)\ p d n x) 1 ^ p . The information entropy 
of ip (or p) is defined as 

S{$) = -J \og(p(x))p(x)d n x. (1) 

It measures the localization of the state in configuration space. A high entropy implies a 
low spatial localization and vice versa. Likewise, one can consider the wave function in 
momentum space, defined by the Fourier transform of i/j, that is 

{T^){x) = Je 2 ^(y)d n y (2) 

for if) integrable. (The normalization of T corresponds to using units 2nh = 1.) We will 
often use the notation ip for the Fourier transform of ip. Again, its information entropy S(ip) 
is a measure of its momentum space localization. 

As shown by Hirschman |]J in one dimension and by Bialynicki-Birula and Mycielski 
in the n-dimensional case [Q, the bàsic uncertainty relations of position and momentum in 
quantum mechanics can be derived from the following sharp bound in L 2 (R n ): 

S(i(>)+S(j>)>n(l-log2). (3) 

Indeed, this inequality puts a bound on the maximum localization in phase space and, in par- 
ticular, it can be shown to imply the uncertainty relations of Heisenberg (Weyl-Heinsenberg 
inequality) [0,0. As stressed by Deutsch 0, entropic uncertainty relations among observ- 
ables are a more faithful expression of the quantum-mechanical uncertainty principle than 
the customary generalized Heisenberg relations. (See also [f§-|6|] for further details and ap- 
plications.) 

The equality in (^) is reached by any Gaussian function and moreover these are the 
unique minimizers |7[]. Since the Gaussian can be taken centered at the origin, the same 
sharp bound holds in the subspace of even functions. Less obvious is the value of the sharp 
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bound in the subspace spanned by the odd functions, i.e., i/j(—x) = —íp(x), as well as the 
form of the associated minimizing functions. Such question would arise, for instance, in the 
case of two electrons in a triplet spin state since the relative coordinate wave function must 
be odd. In this paper we will address this problem in the one dimensional case, n — 1. For 
future reference, we will denote the functional S(ip) + S(ip) by S(ip) and the subspace of the 
odd functions in £ 2 (R) by 7Y_. Thus, we seek to find the infimum of S in the space T~í-, and 
also to establish the form of the possible minimizers, or, more generally, of the minimizing 
sequences. 

Quite likely the problem just raised is non trivial if treated in a fully rigorous math- 
ematical manner. In 1957 it was noted by Hirschman (in the one dimensional case) that 
the l.h.s. of (^) is non negative; this result follows from the classical Hausdorff-Young in- 
equality (see e.g. 0), he then conjectured that the sharp bound was attained by Gaussian 
functions [p]]. It was not until 1975 that Beckner ||, motivated by preliminary results of 
Babenko ||, established the necessary sharp version of Hausdorff-Young inequality from 
which Hirschman-Beckner inequality immediately follows. On the other hand, the prob- 
lem of finding sharp bounds in restricted classes of functions, such as linear subspaces, seems 
to have deserved less or not attention at all. Given the difficulty of the problem, we have 
adopted here an exploratory approach in order to gather "experimental" information on the 
minimizing function, namely, by numerically minimizing the entropy functional. From the 
point of view of rigorous mathematical results, this procedure can only yield upper bounds 
on the sharp bound, nevertheless it can provide useful insights and help to make reasonably 
founded conjectures on the form of the minimizers. Such conjectures are presented below. 

Let us briefly describe the numerical procedure used. We have considered the expan- 
sion of the elements of L 2 (R) in terms of the orthonormal harmònic oscillator basis <j) n (x) = 
h n (x)e~ nx2 , where h n (x) are the associated Hermite polynomials. Thus ip(x) = Y^=o a n4>ri{x) 
in the mean. In this basis the Fourier transform takes the simple form ip{x) = J2%Lo i n o, n 4> n (x). 
The entropy functional S is then transformed into a function of the complex coefficients a n 
and the problem consists in minimizing this function with respect to Ü2n+i, n = 0, 1, 2, . . ., 
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keeping a2n = and Yl^Lo \ a n\ 2 finite. To address this problem we actually consider the 
following iV-dimensional subspace of 7i_ 

AT-1 

4>( x ) = H a 2n +i4>2n+i(x) , (4) 

n=0 

for iV as large as possible, then make use of Standard numerical algorithms to look for the 
minimum of S(ip) in this space. The numerical minimization algorithms become less efficient 
as iV increases, thus implying a maximum admissible value for N in practice. The largest 
space used was that corresponding to N = 128, which, of course, yielded the best (i.e., 
the lowest) entropy, namely, S(ip) = 0.61370581. This number, as well as the minimizing 
functions itself, is only very weakly dependent on the minimization method used (e.g. a 
steepest descend or a simplex algorithm), the precise value of N and the initial conditions 
used. Also, we have checked that the Gaussian minimum 1 — log2 is correctly reproduced 
if even as well as odd functions are allowed. It turned out that imposing the conditions 
ip* = ip and x/j = did not result in an increase of the entropy. Analogous restrictions can 
be imposed on the Gaussian minimizer in the subspace of even functions. The minimizing 
function (for iV = 128 and the above mentioned restrictions) is shown in Figure 1. 

Motivated by the numerical results, we define the following two one-parameter families 
of functions, 

$ a (ar) = J2(-l) n e- na2x2 e- n{x ~ n -^ 2/a \ 

& a (x) = J2 (-l) n e- na2{n+ ^ 2 e- 7T{x ' n -^ 2/a2 , (5) 
neZ 

where the parameter a takes positive vàlues. Note that $ a and <&' a are two unrelated func- 
tions; the symbol ' is used to distinguish them and it does not denote a derivative. Our 
preliminary ansatz is that the small a limit of $ a (or equivalently of corresponds to a 
minimizer of S in 7i_. Under this assumption, the numerical curve shown in Figure 1 would 
be a regularized approximation to the small a limit of $ a . In fact, the numerical curve 
coincides almost perfectly with <3> a or <è' a for a = 0.29. The parameter a plays the role of a 
regulator in eqs. @, similar to value of iV in eq. (|). 
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For convenience we will refer to ($ a ) and (<&„) as sequences since it is always possible to 
choose a positive sequence (a n ) with lim^oo a n = so that ($ a „) is a sequence in the usual 
sense. Strictly speaking the limits as a — > of the (<& ) or ($[J do not take place within 
L 2 (R), i.e. in norm. Indeed, their point-wise limit is except at the points x n — n + |, 
where they take the value 1, whereas their norms (H^alb) and (H^'J^) converge to l/y/2, as 
will be shown below. On the other hand the limit of (||$ a — 3^1 1 2) is 0, thus both sequences 
(<& ) and ($„) become equivalent for small a. 

We will introduce the following notation. Let V be a normed vector space, and let (x a ) 
and (y a ) be two sequences in V (in the sense a — > and a taking positive vàlues). We 
will say that they strongly approach each other if lim a ^ ||x — y a \\ = 0, and this will be 
denoted bj x = y oi x a = y a . Let us remark that the sequences are not assumed to be 
Cauchy sequences, hence nothing is implied for the limits of ||x 0l — x 02 || or ||x ai — y a2 \\ as 
ai and a 2 independently approach 0. From the triangle inequality it follows that this is an 
equivalence relation. Furthermore, if x a = y a , it follows that lim a _>o(| \x a \ \ — \ \y a \\) = 0, since 
Kigali — Hï/olDI < \\ x a — UaW- With this notation $ a = &' a in L 2 (R). This is proved in 
Lemma 1 below. 

The word "limit" applied to the sequences (<& ) and ($„) is used here only in an improper 
sense. The strict statement, if correct, would be that (<I> ) is a minimizing sequence in 7í-, 
that is, one that approaches the infimum of S in H-. A calculation, to be discussed in more 
detail later, shows that lim a ^ «S($ a ) = 2(1 — log2), therefore we formulate the following 
conjecture: 

Conjecture 1. The infimum of the functional S in 7Y„ is 2(1 — log2). 

Our best numerical value for S (0.61370581) is only slightly above 2(l-log2) (0.61370564). 
Let us make some remarks on the form of the assumed minimizing sequence and its improper 
limit $o- Both $ a and & a are odd and real functions and moreover jF$ a = i<&' a and vi ce 
versa, thus jF$ a = i$ a in L 2 (R) (and also point-wise). $ is formed by a set of localized 
states arranged antisymmetrically around and distributed equidistantly through the real 
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line. The small scale structures (the so called ultraviolet region in physics) are narrow Gaus- 
sian functions, namely, e ~ 7í ( x - x ^) 2 1 - 2 _ Likewise, the large scale structure (infrared region) is 

2 2 

a wide Gaussian function centered at the origin, i.e, e~ wa x . We will refer to this overall 
arrangement as a singular "bi-Gaussian" function. As it will be shown below, the double 
Gaussian structure of the minimizer <£> is directly responsible for the fact that 2(1 — log2) 
is twice the infimum of S in L 2 (R), which is saturated by a (simple) Gaussian function. As 
we will recali in a moment, minimizing S is equivalent to maximize the Fourier transform 
operator. For any linear operator from L p (R n ) into L q (R n ) with a Gaussian kernel, Lieb has 
shown that the (unrestricted) maximizer, if any, must be a Gaussian function 0. The sub- 
space 7Y_ does not contain such functions, yet it seems that nevertheless a kind of Gaussian, 
to wit, a bi-Gaussian function, is the maximizer also in this case. It is remarkable how even 
the seemingly simple constraint ip £ yields a quite unexpected result, which however is 
full of structure. 

The sharp bound n(l — log2) on the entropy S follows from computing the norm of the 
Fourier transform operator T considered as a linear operator from the space L p (R n ) into 
its dual L 9 (R n ), with p- 1 + q- 1 = 1 and 1 < p < 2 < q 0. Indeed, we can define a new 
functional as 

S q vanishes at q — 2, since T is unitary in L 2 (R"). The functional S q is related to the S by 

S(rl>)-- 



(7) 

q=2 



dq 

where the derivative is a right derivative. Following Hirschman's argument, let K q (V) denote 
the norm of the operator T restricted to a subspace V of L p (R n ), i.e. 

mf{S q (iP), ^jeV} = -\ogK q (V). (8) 

From eq. (|^) and using K 2 (V) = 1, it follows 

H K (V\ 

(9) 



inf{5(^), ^ e V} = -4 dK i} V ï 

dq 



q=2 



As first proved by Beckner ||, the infimum in L p (R n ) is reached by Gaussian functions and 
thm K q (D>(R n )) = (pVPq-Vi)"/ 2 . In view of Conjecture 1, it is natural to make the stronger 
assumption 

Conjecture 2. The norm of the linear operator T from the space 7i p _ of odd functions 
of L P (R) into L q (R), with 1 < p < 2, is K q (Ji, p _) = p 1 /Pq~ 1 /i . Correspondingly, the infimum 
of S q in H p _ is i logç - 1 logp. 

Conjecture 1 follows from this one. A calculation to be detailed below shows that the 
sequences ($ a ) and yield the value S q = - logg — - logp as a goes to 0, thus, according 
to this conjecture, they are minimizing sequences also for S q in TíZ.. 

Conjectures 1 and 2 settle the point (or, more properly, open the question) of the infimum 
of S and S q in Tí 1 ^. As noted, we do not expect a strict minimizer of S q to exist and we have 
instead to consider minimizing sequences, i.e. such that lim^oo S q (if) n ) = inf{<5> g ,in TL P _}. 
To address this point and the related problem of uniqueness, and also to give further support 
to the conjectures, we will now turn to a more detailed study of the bi-Gaussian ansatzs $ a 
and <fe' a and their admissible generalizations. 

Let D denote the class of distributions, do(x) of the form 

d {x) = ^2 b n 5(x - x n ) , x n = x + nr, \b n \ = b , r, b > , (10) 
neZ 

for some Xq, r and b. Here S(x) is Dirac's delta function. Further, let ip2 be a function in 
L P (R 2 ). Then, for each a > we will associate to tp2 two functions ipi and ip[ in L P (R) by 
means of the relations 



0*0 = / Mv)^2 (ax, d y=J2 ^2 Ux, x ^ Xn ^ , (11) 

ngZ 

i'iix) = J do(v)ifa (ay, — -^j dy = ^ b n ^ 2 (ax n , ^" j . (12) 

ngZ 

We will use the notations (d ,?p2)a and (d ,ip2)' a to the denote the defining constructions of 
ipi and if)[ respectively. Rather than state the more general conditions on ip 2 for the above 
definitions to make sense, we will restrict ïp 2 to the Schwartz space of fast decreasing C°°(R 2 ) 



7 



functions, on which the tempered distributions are defined. This space will be denoted by T, 
and will be considered as a subspace of L P (R 2 ). It has several useful properties: it is dense 
in L P (R 2 ), is invariant under Fourier transform, and their elements are sufíiciently regular 
for our purposes, in particular, the defining series of ipi and ip[ exist and are absolutely and 
uniformly convergent for given a. 

Both definitions are related by using either ax or ax n as the first argument of tp2- We 
will be interested throughout in the limit of small positive a. In this limit, and for each n, 
the terms faiax, {x — x n )/a) and ^(axn, (x — x n )/a) vanish unless of order a, thus 

both definitions become equivalent. More precisely, they strongly approach each other as a 
goes to 0, i.e. 



This statement is meaningful since the normalizations of ipi an d ip[ are well-defmed; they 
have finite p-norm as a goes to 0. This is proved in Lemma 1 and Proposition 1 below. 

Perhaps the best way of understanding the constructions ipi and ip[ is to consider the 
case of a separable function ^(x, y) = a(x)/3(y). The definition of ipi corresponds to make a 
convolution of d with (3(x/a) and then multiply by a(ax), whereas in tp[ the multiplication 
is performed in the first place and the convolution is done next. In the limit of small a 
both operations commute, i.e., a{j3 * d ) = j3 * (ad ). The function a describes the large 
scale profile of the function, whereas (3 gives the small scale structure of ipi. For arbitrary 
functions tp2, which can always be considered as a linear combination of separable ones, 
those roles are played by ^{x, 0) and ip 2 (0,x), respectively. 

The bi-Gaussian ansatzs $ a and $^ are of the form (do, ^2)0 and (d ,íp 2 )' a with 



(doM a =(doM' a inLP(R). 



(13) 




(14) 



exp(-7r(x 2 + y 2 )) , 



(15) 



i.e. b n = (— l) n and x n — n + |. 
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The small a limit of fa or fa does not take place in L P (R), nevertheless, after an appro- 
priate renormalization, there is a weak limit as a distribution, namely 

lim —fa(x) = lim —fa-, (x) = Kdn(x) , (16) 

where the constant K = J fa(0, x)dx. This follows from considering the integrals / fa(x)f(x)dx 
or / fa 1 (x)f(x)dx for an arbitrary test function /, after the change of variables x — > ax + y. 

Let us compute the p-norm of fa in the limit a — > 0. Since in this limit the overlap 
among different terms of the defining series of fa is negligible, for each value of x at most 
one term of the series is relevant. This can be formulated as follows. For given d in D and 
each integer n, let /„ be interval [x n — \r,x n + \r) and let (p n (x) denote the characteristic 
function of I n . Then 

$i( x ) = S bn ^ 2 ( ax ' X ^(x) , ftiix) = ^2 (ax n , X ^" j <p n (x) (17) 

ragZ ngZ 

represent the truncated functions, obtained keeping only the relevant n for each x, i.e., such 
that x E /„. 

The functions fa, fa 1 and their truncated versions are all equivalent: 

Lemma 1. Let fa E T, d E D and p > 1, then fa = fa 1 = fa = ^4 m L P (R). 

The proof is given in the Appendix. 

We can now compute the p-norms of fa in the limit of small a. 
Proposition 1. Let fa ET, d E D and p > 1, then 

HmlhMlp = i^H^íllp = ^ \\H\p- (18) 
Proof. Due to the previous lemma, it is enough to compute lim a ^ HV'íllp· 
ll^ill^^E / \fa(ax n ,(x-x n )/a)\ p dx 

ZJ In 

/r /2a 
\fa(ax n ,x)\ p dx. (19) 
-r/2a 

neZi 

Due to Lemma 1, the limits of the integral can be extended to ±00. Next, we can use that 
for a Riemann integrable function /, / f(x)dx = lim^o Z) n6 z hf( x o + hn). Thus, 
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lim|K||£ = - / \My,x)\ P dxdy. (20) 

This proves the proposition. 

Corollary 1. Forip 2 ,(p 2 d E D and p > 1, (d ,ip 2 )a = (d>o,<f>2)a in L P (R) if and 
only ifï/j 2 = 02- 

A straightforward calculation shows that the Fourier transforms of ipi and ip[ are related 
to that of ip2 as 

F{do, ifa) a = (d , Tij 2 )'a , F{do, 1p2)' a = (d , TÍJ 2 )a , (21) 

where T denotes the transposition operator, Tip 2 (x,y) = ip 2 (y,x). 

In order to proceed, we will consider admissible only the distributions d in D whose 
Fourier transform d (x) is again in the class D Q , that is 

d {x) = ^2b n ó(x - x n ) , x n = x + nr, \b n \=b, (22) 
neZ 

for some Xq, f and b. The admissible distributions span the space Dq := Z) njF _1 D . Then, 
recalling that T is a bijection in T , Lemma 1 and Proposition 1 apply to ípi and ■íp' 1 . This 
immediately leads to 

Proposition 2. For ip 2 eT, d G Dq and 1 < p < 2, p' 1 + q^ 1 = l, 

lim<S 9 (V>i) = S q (^ 2 ) + C q , C q = log (t^j , (23) 
limS(^) = S{^ 2 ) + C, C = - log (rf) . (24) 

In both cases, the first term depends only on ip 2 and the second one only on do- Furthermore, 
if ip 2 (x,y) happens to be separable as a(x)/3(y), the entropies also split as the sum of the 
entropies corresponding to the infrared part a plus the ultraviolet part j3. Let us denote 
by 7í(do) the (improper) subspace of L 2 (R) spanned by the functions ipi, for given d , in 
the limit of small a. From these formulae follows that the minimum entropy in 7ï(do) 
corresponds to those ip\ associated to a Gaussian ip 2 , i.e. tpi must be a bi-Gaussian function. 
Therefore, 
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Corollary 2. Under the same assumptions as in Proposition 2, the following bounds are 
sharp in H(d ) and are attained by ip2 Gaussian. 

limS 9 (V>i) > -logg- -\ogp + C q , (25) 

a^O q p 

limS(V>i)>2(l-log2) + C. (26) 

On the other hand, S q (ipi) is bounded from below by its infimum in L P (R), namely, 
i logg — ^ logp, hence it follows that any do in Dq must satisfy the following inequalities 

C>-l + log2. (27) 



C q > -— logg + — logp (Kp<2) 



The distribution d in eq. (|Ï4|) satisfies Tdo = id . This is easily proved using Poisson's 
summation formula S ng z exp(ï27rnx) = S n6 z + n )' Thus this do is admissible and the 
sharp bounds apply with r = r = b = b = l, i.e., with C q = C = 0. This is consistent with 
Conjecture 2 since $ a is a bi-Gaussian odd function. Moreover, $ a = i$' a and <È>' a = i$ a . 
This follows from fafei y) — ^{y, %) for 1P2 i n eq. (|I5). 



As a further check of Conjecture 2, let us show that, for the admissible do, the functional 
5 g is stationary at ipi, when ?/>i is a singular bi-Gaussian function. 
Lemma 2. Leí <i G Dq, then b 2 /r = b 2 /r. 

Proof. This follows from using that T is unitary in L 2 (R n ), thus, for any ip2 G 

/^2 z^2 

-\\ih\\l = limllVilli = lixn II^Ml = -||TÍ 2 ||2 = -H^llt (28) 

Proposition 3. Leí ^2 be a Gaussian function, and d G Dq, then S and S q are 
stationary at ipi = {do, 1/12)0, in the limit a — > 0. 

Proof. Let us consider a first order variation of ip. The first order variation of S q in 
L p (R n ) is easily computed from its definition, yielding 

5S q (ilj) = Re J 5ip*(x) (G p - T~ X G^F) ip(x)d n x , (29) 

where we have introduced the non linear operator G s 
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Since any Gaussian function is a minimizer of S q , it follows that (G p — J r ~ l G q J z )ip 2 vanishes 
identically when ip 2 is Gaussian. On the other hand, since ip\ bi-Gaussian is a minimizer in 
the subspace TC(do), SS q (ipi) will also vanish if the variation 5ip* is in this subspace. What 
has to be proved is that in fact 8S q {ip\) vanishes under arbitrary variations, in the limit of 
small a. This is equivalent to prové that (G p — T^GqT^x vanishes when ipi is a singular 
bi-Gaussian function. From arguments entirely similar to those used to establish Proposition 
1, it follows 

T 

G s (d , Va)a = ^{d ,G s ip 2 ) a , (31) 
from where it is finally obtained 

(G v -T- l G q T){d ü M a = (d , (L Gp -^jr-i Gq j^^ . (32 ) 

Now, from Lemma 2, r/b 2 equals r/b 2 ; this quantity factors out and the right-hand side 
vanishes for ip 2 Gaussian. This completes the proof. 

Let us note that the inequalities (|27|), as well as Lemma 2 are statements on the space 
Dq only, independent of the construction (do,ip 2 ) a . This construction, however, defines a 
regularization of do which has proven useful to establish properties in Dq. 

It is also interesting to note that a similar construction to that of $ a can be carried out 
for the space of even functions, using the same ip 2 given in eq. ( |H)[ ) and the distribution 

d o( x ) = E^-íí), (33) 
neZ 

which satisfies dç, = d . All the previous arguments apply here and the same sharp bound for 
S q in Tí(d ) is obtained as for the case of d odd. Of course, the corresponding bi-Gaussian 
is known not to be a minimizer of the even functions subspace and at most it can be a 
relative minimum. 

An immediate question is that of the uniqueness of the minimizing sequence. To study 
this point we have first to consider the symmetries of S and S q . In L 2 (R n ), the functional S is 
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invariant under: (i) multiplication by a non-vanishing complex constant, ip{x) \— > Xtp(x), (ii) 
afiine regular transformations, ip(x) i— > i])(Ax + 6), (ra) complex conjugation if)(x) i— > i/)*(x), 
and (ró) Fourier transform, ^(x) i— > ■0( ;r )· The functional <S g in L P (R") is invariant under 
the transformations (i)-(iii) above, whereas under Fourier transform it satisfies Sp^Fip) = 
—S q (ip), provided the corresponding norms exist. In 7í p __ translation invariance does not 
exist and linear transformations consist only of dilatations. The improper minimizer $o 
is invariant under complex conjugation and Fourier transform but breaks dilatation and 
normalization invariances (as also does any non trivial function in L 2 (R™)). In physics 
language, these two symmetries are spontaneously broken. A similar statement can be 
made for the Gaussian minimizers in L p (R n ). Let us remark, however, that the group of 
symmetries generated by transformations (i)-(iv) does not act transitively on the set of 
minimizers of S in L 2 (R), since complex (rather than real) affïne transformations would 
needed to connect two arbitrary Gaussian functions. Likewise, the previous symmetries 
applied to $ do not exhaust the set of minimizers, and in fact the symmetry group in 
T~í(d ) is even larger; e.g. two independent dilatations applied to d and ip2 still define an 
symmetry transformation which acts effectively on 7í(d ) (always meaning in the limit of 
small a for which H(d ) has been defined). 

In passing, note that under a dilatation do(x) — > fj, 1 ^ 2 do(fj,x), (/i positive), the quantities 
b, b, r and f scale as /i~ 1//2 6, yU 1//2 &, fi~ 1 r and fir, respectively. Thus the quantities r/b 2 , C q 
and C are dilatation invariant, as they should. 

As noted above, the sequence (<& ) is not convergent in L P (R) and it cannot be made 
convergent by a suitable (a-dependent) renormalization of $ a since its weak limit is the 
singular distribution do- Thus it is not a Cauchy sequence; two elements $ ai and $ a2 need 
not be near each other in the strong topology, even for arbitrarily small vàlues of a\ and a 2 . 
That is, the sequence does not even approach itself in the mean and hence a precise definition 
is needed to state that some other minimizing sequence must approach this one. Besides, 
note that given a minimizing sequence, one can apply independent arbitrary symmetry 
transformations for each value of a and still have a minimizing sequence. This implies that 
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a minimizing sequence needs not approach strongly ($ a ) or more generally ((do,ip2)a) for 
fixed (i.e. a-independent) d and ip 2 - 

The numerical calculation shows that the infimum in 7i_ can be achieved in the subspace 
T = +i, whereas that corresponding to the subspace JF = —i is larger. This suggests that 
Fourier transform invariance is not spontaneously broken, that is, that after an appropriate 
dilatation, the minimizer can be brought to the space T = +i (note that T is not invariant 
under dilatations). This is similar to the problem of minimizing S in L 2 (R); a minimizer (a 
Gaussian) is not necessarily an even function, but it can be brought to one after a suitable 
translation. In the case of TiF-, not existing a true minimizer, it is important to specify in 
which sense the minimizer must satisfy the condition Tip = (assuming our conjecture 
of unbroken Fourier transform invariance to hold). One can expect that the condition is 
satisfied in the weak sense. This is consistent with the fact that a~ 1 ipi weakly converges 
to do- It cannot be expected, however, to hold in strong sense for an arbitrary minimizing 
sequence. This can be seen noting that every centered Gaussian ?/> 2 , together with do in 
eq. (13), would yield a minimizing sequence in TCL. By a centered Gaussian, it is meant 
a function of the form ^(x, = iVexp(— |Ar 2 — \By 2 — Cxy), where N, A, B and C 
are complex numbers, N is non vanishing and the real part of \Ax 2 + \By 2 + Cxy is a 
positive defmite quadratic form. In this case ipi — (do, 1^2)0, is an antisymmetric bi-Gaussian 
function. Then, íp[ — itpi = i(d , Típ 2 ~'4 , 2)ai whose norm (from Corollary 1) does not go to 
unless ip2 = Tfa, and this equality does not hold for an arbitrary centered Gaussian ip 2 - 

Another consideration follows from noting that the information entropy S of do is un- 
defined; each single delta function has entropy minus infmity since they correspond to a 
maximal localization, however, the fact that this localization can occur in any of the points 
x n with equal probability adds a plus infmity to the entropy yielding a undefmed value. 
It follows that the value of S or S q for a sequence in 7i{do) depends not only on its weak 
limit, do, but also on the particular shape of the functions: the true minimizer must be 
Gaussian-like. 
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Perhaps it will be useful to illustrate the situation with an example. Consider the mini- 
mization of the functional F(ip) = F (íp) + Fi(ip) on L 2 (R), where 

F (V) = - \ log((a: - (x)^) 2 )^ , iq(V>) = ((x - (x>^)% , (34) 

and (f(x))j, means / f(x)p(x)dx. F has been adjusted so that it is invariant under di- 
latations, whereas, F\ is minimized by functions as narrow as possible. Therefore, we can 
proceed by classifying the space of functions by their value of F\, and choose the minimizer 
of Fq in each class. A simple calculation, using Lagrange multipliers, shows that the mini- 
mizer is a Gaussian located anywhere and with arbitrary normalization and a well-defmed 
width a. This gives F\ = a 2 /2ir and F = ~ + ilog27r. Next, in order to minimize Fi, 
we should take a — > 0. The infimum of F is then | + |log27r. The absolute minimizer 
does not exist in L 2 (R), but a minimizing sequence must approach in some sense the se- 
quence ip a {x) = exp(— nx 2 /2a 2 ) in the limit a — > 0, modulo normalization and location. 
Furthermore, the corresponding probability density p a {x) must approach the distribution 
S(x), again modulo translations. 

After these considerat ions, we will state our conjecture on the form of the minimizer of 
S q in 7i p _. Essentially, it is that a minimizer must necessarily be a singular bi-Gaussian in 
the space T = +i (in the weak sense) and modulo dilatations. To put this conjecture in 
precise terms, let d denote precisely the distribution in eq. (|Ï4"D and let ip^^x) denote ip(px), 
where /i > and if) G L P (R). 

Conjecture 3. Let (ip a ) be a minimizing sequence for S q in 7i p _ (in the sense of a — ► 
and the parameter a taking positive vàlues). Then, (a) there is a sequence of positiu e numbers 
(p a ) and a sequence of complex numbers (X a ) such that {\ a ip^ a ) converges weakly to do- 
According to (a), let us assume, without loss of generality, that (ip a ) has this property with 
p a = 1 and that it is normalized to unity, \\ip a \\p — 1- Then,(b) there is a sequence of 
centered Gaussian functions (ip2,a) such that the sequence ((do,ip2,a)a) strongly approaches 
(ipa), that is, lim a ^o HV'a - (d ,ll>2,a)a\\p = 0. 
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It is clear that this conjecture is stronger than Conjecture 2. It states that conditions (a) 
and (b) are necessary for a minimizing sequence. On the other hand, assuming Conjecture 2, 
they are sufncient: first note that ip a = <p a in L P (R) guarantees ip a = (p a in L q (R), and thus if 
H^oIIp is normalized to unity and ||^ a ||g has a finite non zero limit, \im a ^ (S q (ip a ) — Sq(<f>a)) = 
0. Further, Proposition 1 was proved assuming ip 2 (x,y) to be independent of a. The danger 
with an a-dependent ip2,a is that, if a 2 (|/ 2 )^ 2ia or a^ 2 {x 2 )^, 2 a do not go to for small a, the 
various terms in the series of ip\ or ipi, respectively, overlap and the proposition does not 
apply. This danger is avoided by condition (a) since ip a is assumed to approach d which 
consists of well separated Dirac deltas. 

In conclusion, we have presented a set of conjectures on the infimum and on the min- 
imizers of the functionals S and S q in the space of odd one-dimensional functions. They 
are based on information obtained through a simple-minded direct approach, namely, a 
numerical minimization. This cannot be made into anything rigorous, since the numerical 
procedure might be lead to a relative minimum, rather than to the absolute one, however 
this possibility seems quite unlikely to us since the numerical result has been checked to be 
stable against details of the calculation, including changes in the initial conditions chosen 
for the minimization. 

Although at first sight the numerical result in Figure 1 seems to be rather irregular, we 
have hopefully shown in this work that in fact it is plenty of structure and regularity. The 
space Tí(do) has proven to have nice properties directly inherited from the map d ® L P (R 2 ) 
into L P (R). The numerical value of the infimum of S in has been understood as an 
approximation to twice the absolute infimum in £ 2 (R) and the numerical minimizer has 
been understood as a (singular) double Gaussian structure. Gaussian functions seem to 
dominate the entropy minimization problem both in the whole space and in the odd functions 
subspace. Likely, these regularities will open the way for a rigorous treatment of the problem 
studied here. 
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APPENDIX A: PROOF OF LEMMA 1. 



Let us begin by proving ip[ = ip[. For p > 1, || || p is a norm, hence, due to the triangle 
inequality 



IK - V'ÍHp < 6 1] (/ 1^2 (a^n, (a: - x n ) / a)\ p dx) 

n,m \J Im / 



(Al) 



By assumption ^ 2 is bounded and fast decreasing at infinity, thus for any s and t positive, 
there is a positive K such that \ip 2 (x,y)\ < K{\ + x 2 ) _s |y| _í . Also, for x £ I m , \x — x n \ > 
r(\n — m\ — |). Therefore, for s > 1/2 and í > 1 



|K-#|| p <&tfa*£(l + a 2 4) 'Ei/ \x-x n \-»dx) 



\ 1/P 



n — m — 



m>l 



neZ 



< 7Ta 



'„í-i 



(A2) 



In the last inequality we have used that the series on n is of order a _1 since (1 + x 2 )~ s 
is Riemann integrable and the x n are equidistantly distributed. The proof of ipi = ipi is 
analogous. 

Since = is an equivalence relation, it remains only to show that ipi = 



Z->/„ 



dx 



< ò p a E / 1^2 (aa; n + a 2 x, x^j - *02 (ax n , x) 



dx . 



(A3) 



Again, dip 2 (y , x) / dy is a fast decreasing function, thus, choosing s > 1/2 and í > 1, 



ip 2 (y + a 2 x,x) - ip 2 (y,x) 



< 



< 



< 



ry+a x 

/ \d z i/j 2 (z,x)\ dz 
Jy 

Jy (1 + X 2 ) Í (1 + ^ 2 ) 



-dz 



a 2 K\x\ 



(l + x 2 Y(l + y 2 Y 



(A4) 
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Ilíi -í'i \\l < WK p a 2p+1 í —^—-àx V t~ ... 

< ifa 2p (A5) 

This completes the proof of the lemma. Note that the conditions imposed on ip 2 are far 
more restrictive than actually needed in the proof. 
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